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We present a detailed study of non-leptonic two-body decays of B mesons based on 
a generalized factorization hypothesis. We discuss the structure of non-factorizable 
corrections and present arguments in favour of a simple phenomenological descrip- 
tion of their effects. To evaluate the relevant transition form factors in the fac- 
torized decay amplitudes, we use information extracted from semileptonic decays 
and incorporate constraints imposed by heavy-quark symmetry. We discuss tests 
of the factorization hypothesis and show how unknown decay constants may be de- 
termined from non-leptonic decays. In particular, we find /u^ = (234 it 25) MeV 
and fo" = (271 ± 33) MeV. 



1 Introduction 

The weak decays of hadrons containing a heavy quark offer the most direct 
way to determine the weak mixing angles of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix and to explore the physics of CP violation. At the same time, 
they are of great help in studying strong-interaction physics related with the 
confinement of quarks and gluons into hadrons. Indeed, both tasks complement 
each other: an understanding of the connection between quark and hadron 
properties is a necessary prerequisite for a precise determination of the CKM 
matrix and CP-violating phases. 

The simplest processes are those involving a minimum number of hadrons, 
i.e. a single hadron in the final state of a semileptonic decay, or two hadrons in 
the final state of a non-leptonic decay. In recent years, much progress has been 
achieved in understanding these processes. Simple bound-state models are able 
to describe, in a semiquantitative jssiay, the current matrix elements occurring in 
semileptonic decay amplitudes cIlII. A factorization prescription for reducing 
the hadronic matrix elements of four-quark operators to products of c urrent 
matrix elements shed light onto the dynamics of non-leptonic processes 113113, 
where even drastic effects had been lacking an explanationjpcfpre. 

More recently, the discovery of heavy-quaxk symmetry Ejil3 and the estab- 
lishment of the heavy-quark effective theory E3l1I have provided a solid theo- 
retical framework to calculate exclusive semileptonic transitions between two 
hadrons containing heavy quarks, such as the decays B D'^*^£i>. Moreover, 
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the heavy-quark expansion has been apphed to—thii. calculation of inclusive 
semileptonic, non-leptonic aacLxare decay rates E3 lj and of the lifetimes of 
charm and bottom hadronscJ~nZl. These developments are discussed in detail 
in the article of one of us (M. Neubert) in this volume. In particular, they lead 
to a precise determination of the CKM matrix element governing the strength 
of 6 ^ c transitions: 

iV^fcl = 0.039 ± 0.002. (1) 

In this article, we shall discuss non-leptonic two-body decays of B mesons. 
The dynamics of non-leptonic decays, in which only hadrons appear in the 
final state, is strongly influenced by the confining colour forces among the 
quarks. Whereas in semileptonic transitions the long-distance QCD effects 
are described by few hadronic form factors parametrizing the hadronic ma- 
trix elements of quark currents, non-leptonic processes are complicated by the 
phenomenon of quark rearrangement due to the exchange of soft and hard glu- 
ons. The theoretical description involves matrix elements of local four-quark 
operators, which are much harder to deal with than current operators. These 
strong-interaction effects prevented for a long time a coherent understanding 
of non-leptonic decays. The A/ = ^ rule in strange particle decays is a promi- 
nent example. Although this selection rule had been known for almost four 
decades, only recently successful theoreiicaJ approaches have been developed 



to explain it in a semiquantitative wayEa~0. The strong colour force between 
two quarks in a colour-antitriplet state has been identified as the dominant 
mechanism responsible for the dramatic enhancement of A/ = ^ processes. 

The discovery of the heavy charm and bottom quarks opened up the pos- 
sibility to study a great variety of new decay channels. By now, there is an 
impressive amount of experimental data available on many exclusive and in- 
clusive decay modes. In many respects, non-leptonic decays of heavy mesons 
are an ideal instrument for exploring the most interesting aspect of QCD, i.e. 
its non-perturbative, long-range character. Since the initial state consists of an 
isolated heavy particle and the weak transition operator exhibits a well-known 
and simple structure, a detailed analysis of decays into particles with different 
spin and flavour quantum numbers provides valuable information about the 
nature of the long-range forces influencing these processes, the same forces 
that determine the internal structure of all hadrons. 

In energetic two-body transitions, hadronization of the decay products 
does not occur until they have traveled some distance away from each otherLj. 
The reason is that once the quarks have grouped into colour-singlet pairs, soft 
gluons are ineffective in rearranging them. The decay amplitudes are then ex- 
pected to factorize into products of hadronic matrix elements of colour-singlet 
quark currents. The factorization approximation has been applied to many 
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two-body decays of B and D mcsonslij'l It relates the complicated non- 

leptonic decay amplitudes to products of meson decay constants and hadronic 
matrix elements of current operators, which are similar to those encountered 
in semileptonic decays. The decay constants are fundamental hadronic param- 
eters providing a measure of the strength of the quark-antiquark attraction 
inside a hadronic state. As some of them are not directly accessible in leptonic 
or electromagnetic processes, their extraction from non-leptonic transitions 
may provide important information. 

In Section ||, we discuss the effective Hamiltonian relevant for decays of 
the bottom quark. In Section ^, we provide a prescription for the calculation of 
non-leptonic decay amplitudes in the factorization approximation. We discuss 
the problems connected with factorization, such as the choice of a suitable 
factorization scale and its possible dependence on the energy released in a 
decay process. In Section |[ we then discuss in more detail the structure of non- 
factorizable corrections. Using the l/N^ expansion, we argue that in energetic 
two-body decays of B mesons a generalized factorization prescription holds, 
which involves two parameters al^ ~ ci+C,C2 and « C2+C,ci depending on a 
hadronic parameter C = 0{1/Nc)- We estimate that the process dependence of 
C is very mild, suppressed as /^.E/mb, where Ai? is the difference in the energy 
release in different two-body decay channels. Some basic aspects of final-state 
interactions are briefly covered in Section ^j. The remainder of this article is 
devoted to a phenomenological description of non-leptonic two-body decays 
of B mesons. In the factorization approximation, the ingredients for such an 
analysis are meson decay constants and transition form factors. In Section ^ 
we collect the current experimental information on meson decay constants. In 
Section we describe two simple models that provide a global description 
of heavy-to- heavy and heavy-to- light weak decay form fartnrs embedding the 
known constraints imposed by heavy-quark symmetry E3ol£j. With only a 
few parameters, these models reproduce within errors the known properties of 
decay form factors and predict those form factors which are yet unknown. Our 
approach is meant to provide global predictions for a large set of non-leptonic 
decay amplitudes. We do not attempt here a more field-theoretical calculation 
of weak amplitudes, nor do we perform a dedicated investigation of individual 
decay channels^ In Sections || and |[ we compare our predictions with the 
experimental data on two-body non-leptonic decays of B mesons, present tests 
of the factorization assumption, and extract the decay constants of and D* 
mesons. In Section |lO|, we briefly address the interesting possibility of i3-meson 
decays into baryons. A summary and conclusions are given in Section |ll|. 

"Attempts of a more rigorous calculation of some particular decay modes are discussed 
in the article by R. Riickl in this volume. 
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2 Effective Hamiltonian 



At tree level, non-leptonic weak decays are described in the Standard Model by 
a single VF-exchange diagram. Strong interactions affect this simple picture in 
a two-fold way. Hard-gluon corrections can be accounted for by perturbative 
methods and renormalization-group techniques. They give rise to new effective 
weak vertices. Long-distance confinement forces are responsible for the binding 
of quarks inside the asymptotic hadron states. The basic tool in the calculation 
of non-lcptonic amplitudes is to separate the two regimes by means of the 
operator product expansion (0PE)E3, incorporating all long-range QCD effects 
in the hadronic matrix elements of local four-quark operators. This treatment 
appears well justified due to the vastly different time and energy scales involved 
in the weak decay and in the subsequent formation of the final hadrons. 

Integrating out the heavy W-boson and .top-jquark fields, one derives the 
effective Hamiltonian for b ^ c,u transitionsE3'L£l: 

+ Kf, [cl{^i) Qf + C2{^l) Qf] + h.c.j 

+ penguin operators. (2) 

It consists of products of local four-quark operators renormalized at the scale 
/i and scale-dependent Wilson coefficients Ci(/z). Vcb and Vub are elements of 
the quark mixing matrix. The operators Qi and Q2, written as products of 
colour-singlet currents, are given by 

Qf = [{d'u)v^A + {s'c)v-a] {cb)v-A , 

Qf = {cu)v~A {d'b)v^A + (cc)y_A {s'b)v-A , 

Qf = [{d'u)v-A + {s'c)v^a] {ub)v-A , 

Qf = iuu)v-A {d'b)v-A + {uc)v-A {s'h)v-A , (3) 

where d' and s' denote weak eigenstates of the down and strange quarks, re- 
spectively, and {cb)v-A = C7^(l — 75)6 etc. Without strong-interaction effects 
we would have ci = 1 and C2 = 0. This simple result is modified, however, by 
gluon exchange: the original weak vertices get renormalized, and new types of 
interactions (such as the operators Q2)y^e induced. Not explicitly shown in 
are the so-called penguin operators e3. Since their Wilson coefficients are 
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very small, the corresponding contributions to weak decay amplitudes only be- 
come relevant in rare decays, where the tree-level contribution is either strongly 
CKM-suppressed, as in TT, or where matrix elements of the Qi and 

Q2 operators do not contribute at all, as in B — > K*'y and K'~'ip. 

The Wilson coefficients Ci(/i) take into account the short-distance cor- 
rections arising from the exchange of gluons with virtualities between niw 
and some hadronic scale /i, chosen large enough for perturbation theory to 
be applicable. The coefficients ci(/i) and C2(/i), for example, arise from the 
hard-gluon exchanges shown in Figure 0. The effects of soft gluons (with 
virtualities below the scale ^) remain in the hadronic matrix elements of 
the local four-quark operators Qi. In the course of the evolution from mw 
down to the scale there arise large logarithms of the type [as \ii(mw / , 
which must be summed to all orders in perturbation theory. This is achieved 
by means of the renormalization-group equation (RGB). The combinations 
c±(m) — ci(/^)±C2(/i) of the Wilson coefficients have a multiplicative evolution 
under change of the rcnormalization scale. They can be obtained from the 
solution of the RGE 

(M|--r±^c±(/i) = o, (4) 

with the initial condition c±(rnvF) = Ij corresponding to ci{mw) = 1 and 
C2{mw) — 0. The quantities T± in are the anomalous-dimensions of the 
operators {Qi ± Q2)- At one-loop order, they are given byOH 

r. = 7±^; 7. = 6(±i-i-), (5) 

where Nc = 3 is the number of colours. To leading logarithmic order (LO), 
the solution of the RGE is 

''^^^^=[^M^) ' Po^^N^-^Uf, (6) 
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where Pq is the first coefficient of the (3 function, and n f the number of active 
flavours (in the region between mw and /i) . The physical origin of the enhance- 
ment of c_ {fi) can be traced back to the attractive force between two quarks in 
the colour- antitriplet channel of the scattering process b + u ^ c + d. Similarly, 
in the colour-sextet channel the force is repulsive, leading to c+(/i) < 1. 

The leading logarithmic approximation can be improved by including the 
next-to-leading (NLO) corrections of order as[as\n{mw / ^J)\^ ■ The result is 

where the coefficients R± are given byi0 is the two-loop coefficient of 
the (3 function) 
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In Table we show the values of ci(mb) and C2{mh) obtained at leading 
and next-to-leading order. The evolution of the running coupling constant is 
done at two-loop order, using the normalization as{mz) = 0.118 ± 0.003. For 
comparison, we show in Table || the values obtained at a lower scale, which 
may be relevant to charm decays. The negative sign of C2(/i) has important 
consequences for decays in which there is significant interference between ci 
and C2 amplitudes. 

In a typical hadronic decay process such as B'^ q~ , there are several 

mass scales involved: the hadron masses, the quark masses, the energy release, 
etc. Thus, there is an uncertainty in the choice of the "characteristic scale" of 
a process. In principle this is not a problem, since the products of the Wilson 
coefficients with the hadronic matrix elements of the local four-quark operators 
are scale independent. In practice, however, one often employs simple model 
estimates of the matrix elements, which usually do not yield an explicit scale 
dependence that could compensate for that of the Wilson coefficients. Instead, 
these model calculations are assumed to be valid on a particular scale. A re- 
lated technical problem is that whereas the operator evolution from mw down 
to m}j can be calculated in a straightforward way, a more complicated scaling 
behaviour is expected below the mass of the h quark, where all the other mass 
scales start to become relevant. There have been attempts to account for the 
scaling in the region < to;, by summing logarithms of the type [as ln(mb//Lt)]" 
and [cks ln(i?//i)]", with E — 0{mi,) being the larae-energy of a light parti- 
cle produced in a two-body decay of a _B meson OO. However, at present 
the treatment of such corrections is still associated with large uncertainties. 



6 



Table 1: Values of the Wilson coefficients at the scale nii, = 4.8 GeV, both at leading (LO) 

and next-to-leading (NLO) order 



as{mz) 






0.115 
0.118 
0.121 


1.102 -0.239 
1.108 -0.249 
1.113 -0.260 


1.124 -0.273 
1.132 -0.286 
1.140 -0.301 



Table 2: Values of the Wilson coefficients at the scale nic = 1.4 GeV 



as{mz) 


cLO(m,) cLO(m,) 




0.115 
0.118 
0.121 


1.240 -0.476 
1.263 -0.513 
1.292 -0.556 


1.313 -0.576 
1.351 -0.631 
1.397 -0.696 



We shall thus stay with the conventional choice /i = 0(mb) adopted in most 
previous work on non-leptonic B decays. 

3 Factorization 

In weak interactions, a meson (or meson resonance) can be directly generated 
by a quark current carrying the appropriate parity and flavour quantum num- 
bers. The corresponding contribution to ^.dccay amplitude factorizes into the 
product of two current matrix elements 0L3. As an example, consider the 
transition 7r~. The factorizable part of the amplitude is given by 

Afact = Kb «i (^"1 idu)A I 0) {D+\ icb)v m . (9) 

The coefficient ai will be discussed below. The transition matrix 

element is of the same type as that encountered in the semileptonic decay B^ —f 
D~^£~i'. It can be determined using data on semileptonic-decays together with 
theoretical arguments based on heavy-quark symmetryll3. The amplitude for 
creating a pion from the vacuum via the axial current is parametrized by the 
decay constant /tt and is proportional to the momentum of the pion: 

(7r"(p)|J7^75M|0) = i/^p^. (10) 
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Thus it seems natural to assume that the amphtude for energetic weak decays, 
in which the directly generated meson carries a large momentum, is dominated 
by its factorizable part. This assumption can be substantiated by a more 
detailed analyses of the kinematic situation in the above decay process: a fast 
moving (ud) pair created in a point-like interaction, with both quarks leaving 
the interaction region in the same direction and with a velocity close to the 
speed of light, will hadronize only after a time given by its 7 factor times a 
typical hadronization time Thad ~ 1 fm/c. In the above example, this means 
that hadronization occurs about 20 fm away from the remaining quarks. Inside 
the interaction region, the (ud) pair behaves like a colourless and almost point- 
like particle. It does little interact with the-jremaining quarks. Because of this 
intuitive "colour transparency argument" e3, one expects that the factorizable 
part (||) of the decay amplitude does indeed give the dominant contribution 
to the full amplitude. A more formal investigation of this situation, using 
an effective theory for heavy quarks and fast-moving light quarks, has been 
presented in Ref. 72. 

In the following applications, we shall adopt the factorization ansatz also 
for processes in which the 7 factors of the outgoing particles are not necessarily 
large. Examples of such decays are B KJ/ip and B ~* D Dg. In these 
cases, the kinematic argument given above does no longer apply. Nevertheless, 
in a two-body process, the concentration of energy into colour-singlet states 
together with the fact that (axial) vector current matrix elements increase 
with the particle momentum favours the direct current-induced production 
of mesons over more complicated production mechanisms. Only comparison 
with experiment can tell whether factorization is a useful concept also for 
those processes. We will analyse the structure of non-factorizable corrections 
in Section ^, and discuss some tests of the factorization approximation in 
Section ||. 

By factorizing matrix elements of the four-quark operators containedJii 
the effective Hamiltonian (||), one can distinguish three classes of decays EjiIj. 
The first class contains those decays in which only a charged meson can be 
generated directly from a colour-singlet current, as in D+7r~. For these 

processes, the relevant QCD coefficient is given by the combination 

ai = ci(^/) H- Cc2(Ai/) I (class I) (11) 

where C = 1 l^c (Nc being the number of quark colours) , and ^J■f = 0(mh) is the 
scale at which factorization is assumed to be relevant. The term proportional 
to C2(/x/) arises from the Fierz reordering of Q2 operators and factorization of 
the product of colour-singlet currents contained in it. At this point, the colour- 
octet term resulting from the Fierz reordering is simply being discarded. In 
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order to remove this deficiency, one should treat C as a free parameter. This 
will be discussed in more detail in the next section. 

A second class of transitions consists of those decays in which the meson 
generated directly from the current is neutral, like the J/tp particle in the 
decay B — > K J/tp. The corresponding decay amplitude, 

Afact = ^ Vcb K*a a2 {J/^\ {cc)v \0){k\ {sb)v \B) , (12) 

is proportional to the QCD coefficient 

a2 = C2{fJ.f) + C ci{^f) . (class II) (13) 

Due to the different sign and magnitude of the Wilson coefficients (see Table ^ , 
the combination 02 is particularly sensitive to the value of the factorization 
scale and to any additional long-distance contributions (i.e., to the precise 
value of (■)■ The QCD coefficient ai, on the other hand, can be estimated 
quite reliably. 

The third class of transitions covers decays in which the ai and 02 ampli- 
tudes interfere, such as in D^-k^ . Their final state contains a charged as 
well as a neutral meson, both of which can be generated from a current of one 
of the operators of the effective Hamiltonian. The corresponding amplitudes 
involve a combination 

ai -f a; a2 , (class III) (14) 

where x = 1 in the formal limit of a flavour symmetry for the final-state 
mesons, as it is realized in the strongly CKM-suppressed decay B~ tt^tt". 

In principle, there is another type of factorizable contribution to weak de- 
cay amplitudes, which is however significantly differentpirpm the ones covered 
so far: the so-called "weak annihilation contribution" E3£3, in which the de- 
caying heavy meson is annihilated by a current of one of the operators of the 
effective Hamiltonian. For a charged meson, this contribution is proportional 
to oi, while it is proportional to 02 for a neutral meson. In the latter case, 
the weak annihilation is in fact the exchange oi a W boson between the two 
constituent quarks. In a weak annihilation process, the second current in the 
four-quark operator produces all the recoiling final-state particles out of the 
vacuum, which implies a sizable form factor suppression. Therefore, annihila- 
tion amplitudes are expected to be small, and it is commonly assumed that 
they may be neglected for all except some rare processes. This approxima- 
tion is not an essential part of the factorization scheme, however. It would be 
straightforward to include the annihilation contributions if the relevant form 
factors at large time-like values of were known. 
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Let us come back to our first example of a class I decay, collecting all 
the different pieces of the factorizable contribution to the decay amplitude. 
Inserting (|l^) and a suitable form factor decomposition of the hadronic matrix 
element (see Appendix) into (||), we obtain 

Afact(B° ^ D+TT-) = ^ Kb KTd ai h « - ml) F,^'-''{ml) . (15) 

Note that it is the longitudinal form factor Fq which enters this expression, 
whereas (in the limit of vanishing lepton mass) the semileptonic decay ampli- 
tude for 5" D~^£~i' involves the transverse form factor Fi. Fortunately, 
heavy-quarkjSKmmetry provides relations between all the B I?^*) transition 
form factors E^Ej. They allow us to relate Fq to form factors which are easily 
accessible experimentally| 

Note that the hadronic matrix elements of the vector and axial currents 
resulting from the factorization of the matrix elements of four-quark operators 
do not show any scale dependence that could compensate the scale dependence 
of the Wilson coefficients. Strictly speaking, therefore, factorization cannot be 
correct. What one may hope for is that it provides a useful approximation 
if the Wilson coefficients (or equivalently the QCD coefficients ai and 02) are 
evaluated at a suitable scale fif, the factorization point. Because the factorized 
hadronic matrix elements can only account for the interaction between quarks 
remaining together in the same hadron, the Wilson coefficients in effective 
Hamiltonian should contain those gluon effects which redistribute the quarks. 
Thus, we have to evolve these coefficients down to a scale where the gluons are 
no longer effective in changing the particle momenta in a significant way. For 
very energetic decays, the "colour transparency argument" shows that gluons 
with virtualities much below the mass of the decaying particle are ineffective 
in rearranging the final-state quarks LnO, and fJ,f — 0{mi,) seems a reasonable 
choice for the factorization scale. However, for processes in which the decay 
products carry only little kinetic energy, gluons with virtualities well below 
TOfc can still lead to a redistribution of the quarks before hadronization sets 
in. This fact may be used as an argument in favour of a lower factorization 
scale in such processes. On a qualitative level, the connection between the 
factorization scale and the energy release in the final state can be seen from 
Figure ^, where we show the ratio 02/01 as a function of Q;s(/i/). As we will see 
in Section ^, the value preferred by i? Dtt decays is positive and corresponds 
to a rather small coupling, indicating fif — 0{mb). On the other hand, D 

''In addition, in the present case the smallness of the pion mass allows us to use the 
kinematic constraint -Fo(O) = Fi{0) (see Appendix) to relate _Fo to a measurable semileptonic 
form factor. 
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Figure 2: The ratio a2/ai as a function of the running couphng constant evaluated at the 
factorization scale. The bands indicate the phenomenological values of 02/111 extracted from 

B — > Dtt and D — > K-k decays. 



decays indicate a negative value of 02 /ai, corresponding to a significantly lower 
value of the factorization scale. This is in accordance with the fact that in these 
processes the energy released to the final-state particles is only about 1 GeV. 
Moreover, it reflects that for low-energetic processes hadronic weak decays 
exhibit a strong A/ = ^ enhancement, which is significant in D decays and 
even spectacular in K decays. A strict AJ = i selection rule would correspond 
to 02/01 = —1. Because of these intuitive arguments, we are led to abandon 
the "naive" factorization prescription adopted in most previous work, where 
a fixed factorization scale — mi, was taken for all B decays, in favour of a 
more flexible factorization scheme, in which /i/ - or equivalently C in (0) and 
( p!3| ) - is treated as a process-dependent parameter. 

We end this section with a remark on very low-energetic processes such as 
K decays. Factorization of hadronic matrix elements of four-quark operators 
into two matrix elements of colour-singlet currents implies that only those non- 
perturbative forces that act between quarks and antiquarks are taken into ac- 
count. The remaining interactions including, in particular, the gluon exchange 
between two quarks or two antiquarks are treated perturbatively. However, in 
the decays of strange particles the long-distance attraction between two quarks 
in a colour-antitriplet state gives the dominant contribution to A/ = ^ decay 
amplitudes. A detailed analysis of these "diquark effects" can be found in 
Ref. 51. As will be briefly discussed in Section 10, we expect similar effects 
to be important in B decays into baryon-antibaryon pairs. For energetic B 
decays into two mesons, on the other hand, the two quarks (or antiquarks) do 
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not end up in the same final-state hadron, and a pcrturbative treatment of 
their mutual interaction seems justified. 



4 Corrections to Factorization and the 1/Nc Expansion 

Let us now illustrate in more detail the structure of the corrections to the 
factorization approximation (for a similar discussion, see Refs. 59, 60). For a 
given class I or class II two-body decay channel, the effective Hamiltonian (||) 
can be rewritten using Fierz identities in such a way that the quarks are paired 
according to the flavour quantum numbers of the final-state hadrons.*! This 
introduces products of two colour-singlet or two colour-octet current operators. 
The hadronic matrix elements of the latter ones are being neglected in the naive 
factorization approximation. 

Consider again the example of the decay 5° D+tt^. In this case, the 
appropriate form to write the effective Hamiltonian is 

iJeff = ^ I (^Ci(^) + ^) {du)v^A {cb)v^A 

+ 2^^dtaU)v^A{ctab)v-A} + --- , (16) 



where {dtau)v-A = ^7/^(1 — 75)^0^ etc., and we have shown explicitly the 
dependence on the number of colours. For a given decay process, let us now 
define two hadronic parameters as follows: 

(BD,.). , _ {TT-D+\{du)v^Aicb)v-A\B'>) _ ^ 
^1 IMJ — 



7r-\{du)v-A\0){D+\{cb)v-Am 

(BD,tt), . ^ {■K^D+\ {dtaU)v-A {ctab)v-A \B°) 

~ 2{7r-\{du)v^A\0){D+\{cb)v-Am- ^ ' 

They parametrize the non-factorizable contributions to the hadronic matrix 
elements and are process dependent. Here the subscript refers to the colour 
structure of the operator, whereas the superscript indicates the particles in- 
volved in the process. The final-state particle with flavour quantum numbers 
of one of the currents is given last. This superscript is written only when neces- 
sary to avoid confusion. In terms of these new parameters the decay amplitude 
resumes the form given in (^), however, with ai replaced by the new coefficient 



-off 



C2(m) 



ci(A*) + ^) l+s[''''^''\^^) +c,{^,)e^i'''^^\^^). (18) 



'^In all our applications, class III amplitudes arc related to combinations of class I and II 
amplitudes by isospin relations. 
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Similarly, in the class II transition B'^ D'^tt^ we encounter the effective 
coefficient 



off 
0-2 



C2(m) 



Nr 



1 



(19) 



Since we have introduced these parameters without loss of generality, the 
effective coefficients take into account all contributions to the matrix ele- 
ments and are thus fj, independent. In other words, the hadronic parameters 
Si{fi) restore the correct /i dependence of the matrix elements, which is lost in 
the naive factorization approximation. Using the RGE for the coefficient func- 
tions, it is then straightforward to show that the fi dependence of the hadronic 
parameters is, in the leading logarithmic order, given by 





[(-^) 


1 -1- £l(^o) 


+ £si^J■o) 


K+ 


1 

+ 2 


[0-^) 


1 + £i(mo) 


- £8(a*o) 





e8(M) = t; 



1 



1 - — ) £8(^0) + ( 1 - ]^ 



£8(/^o) 



1 

TV? 



;i(mo) 
:i(mo) 



(20) 



where k± = [as{^l) / as{p-o)X'^ ^"^^^ ^ ^'^'^ Mo is an arbitrary normalization point. 
It is interesting to evaluate these expressions assuming that there exists a scale 
Mo = M/ where factorization holds, i.e. where £i(/i/) — 0. It then follows that 



ei(M) 
£8(m) 



1 



-( - — 

2 iv; 



- 1 = 0{l/Nl) 



{K+-n.) = Oil/Nc). 



and expanding in powers of as we obtain 
ei{n)=0{al), £aif^) = 



-— In— + 0(aj. 



(21) 



(22) 



To deduce the dependence on the number of colours, we have used the ex- 
pressions for the anomalous dimension given in (H). Note that the large- A'^c 
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counting rules, ei — 0{1/N^) and es = 0(1/A^c), |-af|e-|independent of the 
assumption that factorization holds at the scale fi / II3'L3. Hence, from first 
principles of QCD, we expect that |ei| <C 1, whereas contributions from eg can 
be more sizable. 

Similar counting rules can be derived for the Wilson coefficients Ci{fi); 
however, they can be obscured by the presence of the large logarithm L = 
\T^{mw/^J')■ In general, we have 

ci{fi) = l + 0{L/N^), C2(m) -0(i/^c). (23) 

At the scale /i = mt, , it is evident from Table ^ that one can consistently treat 
L = 0(1), and thus ci = 1 + 0{1/N^) and Ca = 0(l/iVc). For much lower 
scales, such as = rric, however. Table |^ shows that it is appropriate to take 
L/Nc = 0(1), and therefore a = 1 + 0{1/Nc) and C2 = 0(1). 

Let us now discuss the relation of this general approach to the conventional 
factorization scheme, focusing first on the case of B decays. Naive factoriza- 
tion corresponds to setting £i(mb) = 0, in which case ai and 02 are universal, 
process-independent coefficients, which are simply linear combinations of the 
Wilson coefficients ci{mh) and C2(m;,). In ( p^ and ( p^ ) we have shown how 
these relations are modified by the presence of non-factorizable corrections. 
From the above discussion it follows that non-factorizable corrections are ex- 
pected to be small in class I transitions. In class II decays, on the other hand, 
the contribution proportional to eg is enhanced by the large value of the ra- 
tio |ci/c2|. Hence, as a general rule, we expect sizable violations of the naive 
factorization approximation in class II decays only. Explicitly, evaluating the 
general expressions in the large- iVc limit, we find 

af = cAp) + C2(a^) + ei'''''^\^i)^ + 0(1/7V|) = 1 + 0{L/Nl) , 
af = C2(/.) + c,{p) (^i- + 4'''^'''^)) + 0{L/Nl) . (24) 

In B decays, with L = 0(1) and neglecting terms of order we can rewrite 

this as 

af«l, af «C2(mb) + Cci(mb), (25) 



with 



C ^ + e^^-'^Hm,) . (26) 



Note that, with the same accuracy, the first relation in (|2^) may be replaced 
by af ~ ci(m&) -I- Cc2('7ib). It is important to stress at this point that the 
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naive choice ai = ci + c%JJSfc and 02 = C2 + ci/Nc does not correspond to any 
consistent limit of QCDl3; in particular, this is not a prediction of the 1/Nc 
expansion. Since the parameter eg is of order l/Nc, the two contributions to 
( are expected to be of similar magnitude, and hence C should be considered 
as an unknown dynamical parameter. In other words, the large-iVc counting 
rules of QCD predict that related class I and class II two-body decays, such 
as 5° D+TT- and B° D^n", or B° D+ g- and B" -> D°g°, can be 
described (up to corrections of order using factorization with a single 

phenomenological parameter ( — 0{1/Nc). Strictly speaking, however, this 
parameter will take different values for different class II decay channels. The 
picture emerging from this discussion is equivalent to the concept of using 
a process-dependent factorization scale fj,f, discussed in the previous section. 



This scale is defined such that SsifJ-f ) = 0. Using ( |22| ) one may then calculate 
the corresponding value of esirnb) and thus the C parameter for each process. 

We can go a step further and estimate the expected size of the process 
dependence of the parameter From Figure ^ we know that for the par- 
ticular case of B — > Dtt transitions factorization holds at a high scale, i.e. 
^(Btt.u)^^^^ _ Q ^ thus C « 1/3. If we now assume that in 

other two-body decays of B mesons the factorization scale is lower because 
there is less energy released to the final-state particles, we may use relation 
( p^ to obtain for the change in C: 

. / X ^aJrub) AE AE , , 

37r mb 1 GeV 

where AE is the difference in the energy release in different decay channels. 
In the two-body B decays of interest to us, AE is always smaller than 1 GeV, 
so that AC, becomes a parameter of order l/miy. The corresponding variations 
of C, are of order a few per cent, which is small compared with the value of Q 
itself ~ 1/3 for B Dn decays). 

To summarize this discussion, we repeat that in energetic two-body decays 
of B mesons the large- A^c counting rules of QCD predict that factorization 
with af^ 1 works well for class I transitions, whereas class II transitions 
can be described by a phenomenological coefhcient « C2{mb) + Cci(m(,) 
with C — 0{1/Nc). Based on an intuitive argument, we have estimated that 
the process dependence of C is very mild, since A(^ ^ AE/mb. Thus, we have 
provided a theoretical basis for the phenomenological treatment of using a 
factorization prescription in which the phenomenological parameter ^ is taken 
to be the same for all energetic two-body decays of B mesons. 

Let us briefly also discuss the case of charm decays. Here the large-iVc 
counting rules are different, because empirically L/Nc = 0{1). Then, neglect- 
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ing terms of order we find instead of ( p5| ) tlie relations 

a'f « ci(mc) + C' C2(toc) , 02**" ~ C2{mc) + Cci(mc) , 



(28) 



with 



1 



XDK.n) 
8 




8 



(29) 




In general, class I and II transitions can no longer be described by the same 
C parameter. However, we may again argue that the expected process depen- 
dence of the C parameters is a mild one. In (|2^) now appears the smaller 
charm-quark mass in the denominator, but on the other hand the energy dif- 
ference AE in the numerator is smaller than in B decays. As a consequence, 
we still have of order a few per cent; in particular, then, we expect C « ^. 
That C in ( p5| ) should be treated as a phenomenological parameter, rather than 
fixed to the value C = 1/3 pr|edicted by naive factorization, was the basis of 
the approach of Bauer et al. t2l (see also Ref. 53). It turns out that setting 
C « provides a successful description of many two-body decays of D mesons. 
To some extend, this phenomenological "rule o£, discarding 1/Nc terms" can 
be understood in the context of QCD sum ruleslp, and using more formal con- 
siderations based on the heavy-quark expansionEJ. From the above argument, 
it follows that if C « for one D decay mode, it is expected to be a small 
parameter also for other channels. 

As a final comment, we note that whereas the concept of introducing 
effective, process-dependent parameters a^^ works in most two-body decays of 
B and D mesons, it has to be modified in decays where the final state consists 
of two vector particles (P VV decays). The polarization of the final-state 
particles in such processes is very sensitive to non-factorizable contributions 
and final-state interactions. The ratios between S-, P- and D-wave amplitudes 
predicted in the factorization approximation are affected since non-factorizable 
contributions to the amplitude will, in general, have a different structure for 



non-factorizable contributions to the decay amplitudes are spin-dependent, 
thus affecting the polarization of the final-state particles. Likewise, final-state 
interactions are different for different partial waves. They may change S- into 
D- waves even without changing the total decay rate. A case of particular 
interest is the polarization of the J/tp particle in the decay B — > K*J/tp. Most 
models predict a longitudinal polarization of around 40% (a model that we 
shall introduce in Section ^ predicts 48%), whereas the experimental world 
average isciJ (78 ± 7)%. A recent CLEO measurement, however, givesl£3 (52 ± 
7 ± 4)%. Although a clear picture has not yet emerged, it is possible that 
the first deviations from factorization predictions in B decays will be seen 



different partial waves 



In other words, the matrix elements describing the 
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in polarization data. Indeed, it has been argued that .spin-dependent non- 
factorizable effects are necessary to understand the dataE3 (see, however, also 
the discussion in Ref. 8). 

5 Final-State Interactions 

In the conventional factorization approximation, the scattering of the final- 
state particles off each other is neglected, and all amplitudes are real (apart 
from an arbitrary common phase). Watson's theorem, however, requires the 
amplitudes to have the same phases as the corresponding scattering ampli- 
tudesEj. Moreover, in two-body decays, the final state must be of low angular 
momentum, and we know that S- and P-wave scattering amplitudes in the 
GeV region are, in general, inelastic. 

Let us first consider a scattering cigcnstate, i.e. an eigenstate of the S 
matrix. The phase of the amplitude for the decay into this state, which is in 
general a linear combination of physical final states, is identical to the cor- 
responding scattering phase. But final-state interactions will, in general, also 
affect the magnitude of the scattering amplitudes. By definition, these interac- 
tions occur in a space-time region where the final-state particles have already 
been formed in their ground states, but are still strongly interacting while re- 
coiling from each other. Accordingly, for a very energetic two-body decay with 
a pion in the final state, the "colour transparency argument" of the previous 
section already excludes the possibility of significant final-state interactions, 
since the light quark-antiquark pair will have left the region of strong inter- 
action long before it hadronizes into a pion. In a weak decay involving small 
recoil energies, however, the product kR (particle momentum times the radius 
of the strong-interaction region) is of the same order as the scattering phase 
6, which is 0(1). Therefore, in such a decay rescattering will indeed change 
also the magnitude of the decay amplitude. In if — > 27r decays, for exam- 
ple, rescattering effects lead to an enhancement of the A/ = i amplitude by 
« 30% and to a reduction of the AI = | amplitude by ~ 10%O. Fortunately, 
in very energetic processes such as exclusive B decays, the phase 6 is negligible 
compared to kR, and we may safely assume that the magnitude of the ampli- 
tude for decays into a scattering eigenstate remains unchanged by final-state 
interactions. As a consequence, the relation between the amplitudes At for 
decays into final states i (which are, in general, not scattering eigenstates) to 
the "bare" amplitudes is 




(30) 



where S denotes the strong-interaction S matrix. 
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The S matrix can redistribute the amphtudes into different channels car- 
rying the same quantum numbers, including those channels which were not 



originally coupled to the weak process. To illustrate the consequences of (30) 
we consider, as an extreme example, a two-channel S matrix with maximal 
absorption (5*11 = 5*22 = 0): 

The square root of this matrix is 



(32) 



For a vanishing "bare" amplitude for decays into the second channel {A2 = 
0), the amplitude Ai for decays into the first channel is reduced by a factor 
l/\/2, leading to a 50% reduction in the corresponding branching ratio. The 
second channel, in spite of not being directly accessible through the weak decay, 
obtains an amplitude of the same magnitude. Although this is an extreme 
example, one should be aware of the fact that different decay channels influence 
each other. In particular, decays with small branching ratios may have been 
modified or even caused by a "spill-over" from stronger modes. 

Obviously, using the factorization approximation we can only attempt to 
calculate the "bare" decay amplitudes However, summing over all decay 
channels with the same conserved quantum numbers, the uncertainties con- 
nected with final-state interactions drop out. The reason is that, because of 
( ^ ) and the unitarity of the S matrix, we have 

(33) 

i i 

i.e. this sum of decay rates remains unaffected by final-state interactions. By 
measuring the branching ratios of several related decay channels, one can ex- 
tract the magrLitudes of the corresponding isospin amplitudes as well as their 
relative phasesE3. The same isospin amplitudes can also be obtained from the 
"bare" decay amplitudes calculated in the factorization approximation. Ne- 
glecting inelastic rescattering, one can then extract the QCD coefficients a'^ 
and 02^ from a comparison of measured and calculated isospin amplitudes. 

It is well-known that final-state interactions are important in D decava; 
they must be included in the comparison of theoretical predictions with dataO. 
For B decays, on the other hand, we expect smaller final-state interaction 
effects. Because of the large energy released to the particles in the final state. 
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the magnitude of the isospin amphtudes should not be modified significantly by 
final-state interactions. Moreover, there do not exist any charm resonances in 
the B meson mass region that could lead to strong final-state interactions. As 
an example, consider the decays S — > Dtt. The four-quark operator {du){cb), 
which is part of the effective weak Hamiltonian relevant for B decays, carries 
isospin quantum numbers 1=1 and I3 = +1 and thus can transform a B^ 
meson into a state with / = ^ or / = |, while the B^ can only decay into 
final states with ^ = | . Accordingly, we define the following (complex) isospin 
amplitudes: 

Ay, = {Dtt; / = II = i= {Dtt; / = || iJ^fr IB") . (34) 

They are related with the physical decay amplitudes through 
A(B° ^ D+n-) = ^3/2 + \ll Ai/2 

A(S°^i?V) = y|A3/2-/^^l/2 

A{B- D^TT-) = V3A3/2 . (35) 
Although the phases of the measured decay amplitudes are unknown, we can 



extract from (35) the isospin amplitudes as well as their relative phase: 
= \A{B" D+TT-)\^ + \A{B° D\^)\^ 

-^\A{B- ^D\-)f 

\As/2\' = l\AiB- ^ D\-)\^ 
, , 3|A(ffl-> J+7r-)p-2|Ai/2p-|^3/2p 

Let us assume that the magnitudes of the three physical amplitudes had been 
measured. We could then use these relations to determine the magnitudes of 
the two isospin amplitudes as well as their relative phase. In order to compare 
the results with the amplitudes calculated in the factorization approximation, 
we must assume that the magnitude of the isospin amplitudes remains unaf- 
fected by final-state interactions. In other words, the contribution of inelastic 
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Table 3: Upper limits for the relative phase shifts between the two isospin amplitudes in 
some exclusive hadronie B deeays 



Ratio 


CLEO (90% CL) 1^3/2 - 5i/2\ 


B(B" ^ D07r")/B(B- ^ D^ir-) 
B(B" ^ D°g°)/B{B- D^g-) 
B(B0 D*°7r°)/B{B- D*%-) 
B(bO ^ D*°g°)/B{B- D*°g-) 


< 0.07 < 36° 

< 0.05 < 30° 
<0.13 <53° 

< 0.07 < 36° 



scattering into or from other channels has to be negligible. Then neglecting 
final-state interactions simply amounts to neglecting the relative phase shift 
between the different isospin amplitudes, and we can directly compare the the- 
oretical predictions obtained for |A]^/2| £^nd |j43/2| with the experimental values 
for these quantities. 

Unfortunately, in B decays complete measurements of the class of decay 
amplitudes relevant to the isospin analysis described above are not yet avail- 
able. Still, we can use existing data to obtain upper bounds on the relative 
phase shifts in some interesting cases. In our example, it is the amplitude 
A{B'^ D'^tt'^) that has not yet been measured. However, from the current 
upper limit for the corresponding branching ratio we can derive an upper limit 
for the relative phase shift |(53/2 — ^i/2\- To this end, we represent the second 



relation in (35) as a triangle in the complex plane. An elementary geometrical 
argument shows that the angle between A1/2 and A3/2 is maximal when there 
is a right angle between A1/2 and A{B^ ZJ^tt"). Hence 

• 2,c r N / 9 r(B-) BiB^^D^n^) 

Similar inequalities can be derived for the phase shifts in other decay channels. 
Using then experimental upper limits for, colour-suppressed B decay modes 
obtained by the CLEO CkJlaborationBy, as well as t{B-)/t{B°) = 1.06 ± 
0.04 for the lifetime ratioEIl, we obtain the results shown in Table ||. The next 
generation of CLEO data is likely to produce still smaller upper bounds for the 
various phase shifts, but even the current limits are a strong indication that 
final-state interactions in hadronie two-body decays of B mesons are much less 
important than in the corresponding decays of D mesons. We will therefore 
neglect them in our analysis. 
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6 Decay Constants 



The evaluation of factorized amplitudes requires the knowledge of meson decay 
constants and hadronic form factors of current matrix elements. We shall 
discuss first what is known about decay constants. For a pseudoscalar meson 
P — (9192), we define 

{0\q2l^l5qi\P{p)) ^^fpPf.- (38) 

The decay constant fv of a vector meson P — ((71^2) is defined by'[| 

{0\q27,.qi\Vie)) =e^mvfv- (39) 

The experimental values of the decay constants of the charged pion and 
kaon, as obtained from their leptonic decays l^i^i (7), areE3 

= (130.7 ±0.37) MeV, = (159.8 ± 1.47) MeV . (40) 

For charm mesons, the uncertainties in the experimental values of the leptonic 
decay rates are much larger. For ±hc meson, only an upper bound can be 
deduced from the published dataEj: 

fo < 310 MeV (90% CL) . (41) 

Four groups (WA750, CLEoH BEsEl, and E653E1) have measured the 
H^Vfj, branching ratio, from which the decay constant /d, can be extracted. 
The average value isE3 

fo, (241 ± 37) MeV . (42) 

Later, we shall compare this result with an independent determination of f^^ 
from non-leptonic decays. 

The decay constants of light charged mesons can also be obtained from the 
1-prong hadronic decays of the r lepton. Denoting the relevant CKM-matrix 
element by Vij , and neglecting radiative corrections, we write the corresponding 
decay width as 

r{r-^M-.^) = ^Gl\V.,\'fh (l-^)'(l + ^M^) , (43) 

''For the neutral mesons and ui, we use the current u^piU on the left-hand side of 
( I39I ) and divide the right-hand side by \/2. With this definition and under the assumption 
of isospin symmetry, the decay constants of neutral and charged q mesons have the same 
numerical value. 
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Table 4: Decay constants of light charged mesons extracted from hadronic r decays 



M 


■n K Q K* fli 


Im [McV] 


134 ±1 158 ± 3 208 ±1 214 ± 7 229 ± 10 



Table 5: Decay constants of neutral vector mesons extracted from their electromagnetic 

decays 



V 


Q 


LO 







V'(2S) 


Cv 


1/2 


1/18 


1/9 


4/9 


4/9 


fv [MeV] 


216 ±5 


195 ±3 


237 ±4 


405 ± 14 


282 ± 14 



with bp ~ for pseudoscalar and by — 2 for (axial) vector mesons. We 
use rrir = 177Z..MeV for the mass of the r lepton and = 291.0 ± 1.5 fs 
for its lifetime E3. Moreover, we use the following branching ratios for — > 
M-^^@0: (11.31 ± 0.15)% (tt), (0.69 ± 0.03)% (K), (24.94 ± 0.16)% (g), 
(1.25 ±0.08)% (K*), and (17.65 ± 0.32)% (ai). This leads to the values of the 
decay constants shown in Table ^ Only experimental errors are quoted. The 
results for /^r and fx are in excellent agreement with those given in (^). 

The decay constants of neutral vector mesons can be extracted from their 
electromagnetic decay width, using 

Tiv'^e+e-) = ^—f^cv, (44) 
3 my 

where cy are factors related to the electric chaxge of the quarks that make up 
the vector meson. From the measured widthsEj, we obtain the results shown 
in Table ^. The errors reflect the uncertainty in the experimental data only. 
The value for fg obtained in this way is somewhat larger than that derived 
from r decays. In our analysis below, we take fg — 210 MeV. In all other 
cases, we take the central values shown in the Tables. 

We note for completeness that the decay constants of light neutral pseu- 
doscalar mesons can be extracted from the two-photon deciiv P° 77. Av- 
eragmg the results reported by the TPC/27Ea and CELLOEa Collaborations, 
we find 

= (131 ±6) MeV, = (118 ±5) MeV. (45) 
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Decay constants not yet known experimentally are left as ft'ee parameters 
in our expressions for the branching ratios. However, we will often need an 
estimate for the decay constants of charm mesons, such as D and D*. In 
the heavy-quark limit (nic — > oo), spin symmetry predicts that /u = //j., 
and most theoretical predictions indicate that symmetry-breaking corrections 
enhance the ratio fn* / fo by 10-20%EJ. Hence, we take as our central values 

fo ^ 200 MeV , fo' = 230 MeV , 
fo, = 240 MeV , fo-^ = 275 MeV . (46) 

The value for fo lies in the ball park of most theoretical predictions, whereas 
that for fi)^ corresponds to the central experimental result quoted in (^2|) . 



7 Transition Form Factors 



The most important ingredient of factorized decay amplitudes are the hadronic 
form factors defined in terms of the covariant decomposition of hadronic matrix 
elements of current operators. In particular, we need matrix elements of the 
type 

(M|Q7,,(l-75)6|i?), (47) 

where M is a pseudoscalar or vector meson with mass uim- There are two 
form factors {Fq, Fi) describing the transition to a pseudoscalar particle, and 
four form factors {V, Aq, Ai, A2) for the transition to a vector particle. The 
definitions of these form factors are given in the Appendix. Conventionally, 
they are written as functions of the invariant momentum transfer . To obtain 
reliable theoretical predictions for the form factors is the main obstacle in the 
analysis of hadronic weak decays, once the factorization hypothesis is accepted. 

In the case of the heavy-to- heavy transitions B ^ D and B D* , heavy- 
quark symmetry implies simple relations between the various transition form 
factors. In the heavy-quark limit, they readE3 

F^{q') = V{q^) = A,{ch = A,iq') = !^R±^ ^(^) , 

^o(.^) = ^i(<Z^) = ^^^^e(-), (48) 
mB + niM 2 



where 



ZuiBrnM 
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is the product of the velocities of the two mesons, and the univctsaL (mass- 
independent) function ^{w) is called the Isgur-Wise form factor ES'tS. This 
function is normalized to unity at the kinematic point w = 1, where the two 
mesons have a common rest frame. For the realistic case of finite heavy-quark 
masses, the above relations are modified by corrections that break heavy-quark 
symmetry. They_caii be analysed in a systematic way using the heavy-quark 
effective theory This is discussed in detail in the article by one of us 

(M. Neubert) in this volume. As a consequence, it has become possible to 
extract the B Z?^*^ form factors relevant for all class I transitions considered 
in this article from semileptonic decay data with good precision and in an 
essentially model-independent way. The results are compiled in Ref. 57. The 
use of heavy-quark symmetry constitutes a significant improvement compared 
with earlier estimates of non-leptonic amplitudes, which were based on model 
calculations of the relevant form factors. 

Unfortunately, there has not been similar progress in the calculation of 
hadronic current matrix elements between heavy and light mesons. For these 
transitions, we still must rely on the results obtained using some phenomeno- 
logical model. One such model (below referred to as the NRSX model) wasjfti 
troduced in Ref. 57, where we used the overlap integrals of the BSW modelElEj 
to obtain the form factors at zero momentum transfer, and then proposed a 
specific ansatz for their dependence motivated by pole dominance and the 
relations in (p^). Another modfilQ, which we shall not discuss in detail, has 
been used by Deandrea et al. E3 to perform a global analysis of non-leptonic 
decays comparable to ours. 

In order to get an idea about the unavoidable amount of model dependence, 
which will affect our predictions for class II decays of i? mesons, we will consider 
an alternative to the NRSX model. The main purpose here is not to provide a 
new approach that is more sophisticated; on the contrary, we will make several 
strong assumptions to obtain a model as simple as possible. Yet, we have 
checked that it agrees, within reasonable limits, with the available data on 
semileptonic decays, as well as with most theoretical predictions. Wc thus feel 
confident that this model may also be used advantageously to obtain rough 
estimates for other processes, such as rare B decays and decays of Bs mesons. 



In the case of heavy-to-light transitions, the symmetry relations in (48) 
are no longer valid. We account for this by introducing, for each form factor, 
a function ^{{w) replacing the Isgur-Wise function. There is, however, still a 
connection between various form factors provided by kinematic constraints at 
zero momentum transfer [q'^ = 0), corresponding to 

m?, + mlr 

w = u;„,ax = ^ ■ 50 
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There, the form factors satisfy the relations 



Fi(0) = Fo(0), 

^o(O) = ^^^f±^A,(0)-^^^f^A,(0). (51) 
IniM 2m M 

The normaUzation condition for the Isgur-Wisc function is replaced by the 
inequalities 

U(1)<1, ^A,(1)<1, (52) 

which follow from equal-time commutation relations for a heavy b quark. For 
an estimate of the functions ^i(w) we shall adopt a simple pole model, which 
apart from the masses of some resonances has no tunable parameters. This 
model reproduces the relations (^) in the heavy-quark limit. It will serve for 
an immediate estimate of all form factors in the physical region of q^. As a 
starting point, we write S,i(w) = yj2/{'w + 1) hi{w), where hi{w) will be simple 
monopole form factors. In the heavy-quark limit, the prefactor ^j2/(w + T) can 
be shown to yield an upper bound for the Isgur-Wise functionEj. For heavy- 
to-light form factors, on the other hand, this factor ensures that (,i{w) ^ w~^^^ 
for large w, in accordance with the scaling rules obtained in Ref. 97. In the 



next step we make the strong assumption that the symmetry relations in (48) 
still hold for heavy-to-light form factors, but only at the point = 0. In other 
words, we assume that all the /li(wiiiax) are equal. This is a minimal ansatz 



that preserves the kinematic relations in (51). For the functions hi{w) with 
i = Fi, V, Aq and A2, we take 

fe,(^;) ^ jV ^"^^'-"^' , (53) 

W — Wi 

where TV is a common normalization factor. Note that, according to its defi- 
nition, the value of Wmax in hp-^ is different from that in the other three cases. 
The locations of the poles are at 

= , (54) 

2mBmM 

where Mi denotes the mass of the nearest resonance with the appropriate spin- 
parity quantum numbers. Most of these masses are known experimentally; 
some others (such as the masses of the Be and B* mesons) are taken from 
potential models. To obtain the masses of the 1^ resonances, we simply add 
400 MeV to the masses of the corresponding 1~ states, as suggested by the 
spectroscopy of the light and charm mesons, and use the same values for the 
Ai and A2 form factors. 
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It remains to find an ansatz for hFg{w) and hA^iw). For simplicity, we 
choose to normalize these functions to unity ai w — 1, corresponding to a 
complete overlap of the initial and final states at the same velocity [cf. (|5^)]. 
Moreover, we have to satisfy the condition /i_Fo(?«max) = /lAi(wniax) = 
Thus, we take (for j ~ Fq, Ai) 



^ji'w) = yj _ 1 ' (55) 

1 + r 



where N = 1/(1 + r) is still a free parameter. To fix it, we try to mimic the 
effect of the physical axial vector (1+) pole on the form factor Ai{w) near 
1/7=1. In order to have full correspondence with the pole structure of the 
other form factors in ([5^), the factor ^{w + 1) in (^8|) has to be included. 
Therefore, we require that the derivative of the product ^(w + l)hAi{w) at 
ui = 1 be equal to the derivative of a single-pole form factor with the physical 
pole mass Ma^ ■ The result is 

^_ (niB — mvY 4msmy 




where y is a vector meson. The same value of r is then taken for h (w) . 
With these simple assumptions, all form factors are determined. 

A few tests of this model may be quoted here: The branching ratios for 
the semileptonic decays B D*£p and B D £ D are found ta_be 5.3% 
and 1.8%, respectively. The corresponding experimental values areH (4.64 ± 
0.26)% and (1.8 ± 0.4)%. In the semileptonic decay B q£^^ the values 
for the form factors at = are V = Ai = A2 = 0.26, in accordance with 
the predictions V = 0.35 ± 0J17, Ai = 0.27 ± 0.05 and A2 = 0.28 ± 0.05 
obtained using QCD sum ruleslHa. The expected branching ratio for-this decay 
is 25.5|Kibp, which when combined with the experimental result l£3 B{B ^ 
qIv) = X 10"'' yields |Kb| = (3.1 ± 0.5) x lO'^. The calculated 

p/tt ratio in exclusive_semileptonic B decays is 2.2, as compared with the 
experimental value of l2I 1.4 ± 0.6 

form factor Fi{q^) at threshold with the J5*-pole structure at this point 



i - 




100 


101 



gives gsB'Tv — 0.54 x (0.2GeV//B.) for the BB*tt couphng constant, which is 
of the expected order. 
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8 Theoretical Predictions, Comparison with Experiment and De- 
termination of a^^ and a^^ 

We first discuss the results obtained usina; the NRSX modelB. To obtain pre- 
dictions for the branching ratios from the factorized decay ampUtudes, we 
t{B°) = 1.55±0.04ps andT(B-) = 1.65±0.04 ps for the B-meson hfetimesi 
and Vcb = 0.039 ± 0.002 from (|^). The smah errors in these quantities, as well 
as uncertainties in the decay constants of light mesons, are neglected. Our 
predictions for the branching ratios of the dominant non-leptonic two-body 
decays of B mesons are given in Tables | and ^ The QCD coefBcients af and 
a|*^ (for simplicity called ai and 02 in the tables), as well as the unknown decay 
constants of charm mesons, have been left as parameters in the expressions for 
the branching ratios. For comparison, we show the world average experimental 
results for the branching ratios, as recently compiled in the review article in 
Ref. 102. They are dominated by the CLEO II measurements. The quoted up- 
per limits for the colour-suppressed decay modes also include recent CLEO 
data reported in Ref. 86. In Tables we show for comparison the predic- 

tions obtained using the simple form-factor model described in the previous 
section, including some new decay channels which had not been considered in 
Ref. 57. In the theoretical expression for the class III decay B~ D*'^ai , the 
term proportional to requires the knowledge of the current matrix element 
between the B meson and the pseudovector meson oi. We use the axial cur- 
rent as an interpolating field for the ai particle and employ chiral symmetry to 
relate the corresponding form factors to those of the _B~ — > matrix element. 

Let us now compare in detail the theoretical predictions with the data. 
Due to the uncertainty in the values of the phenomcnological parameters a^^ 
and Oj*^, we first concentrate on ratios of branching fractions in which these 
coefRcients cancel. The comparison of such predictions with data constitutes 
a test of the factorization hypotheses and of the quality of our form factors. 
From the class I transitions listed in Tables || and ||, we obtain the predictions 

B(B° ^ D+TI-) , , 

^(^;-^"/'^,.0.88 [0.89]. (57) 

Here and below we use the NRSX model as our nominal choice, and quote 
results obtained with the new model of Section |^ in parentheses. In the heavy- 
quark limit both ratios become equal to unity, because the spin symmetry 
relates D and D* . Only part of the deviations from this limit is a form factor 
effect; the remainder is of kinematic origin. The experimental results for these 
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Table 6: Branching ratios of non-leptonic B^' decays (in %) obtained in the NRSX modelEJ. 
In the third column, the factors containing not yet known decay const pnta h ave been sup- 
pressed. The last column shows the world average experimental results t^O- Upper limits 
are given at the 90% confidence level 



B° Modes 



NRSX Model 



^cff 
,cff 



1.08 
0.21 



Experimental 
Average 



Class I 



D+TT- 

D+K" 

D+D- 

D+D- 

D+Q- 

D+K*- 

D+D*- 

D+D*- 

D+a^ 

D*+TT- 

D*+K- 

D*+D- 

D*+D- 

D*+g- 

D*+K* 

D*+D* 

D*+D*- 

D*+a7 



0.257 al 

0.020 al 
0.031 al (/d/200)2 
0.879 al (/i3,/240)2 

0.643 al 

0.035 al 
0.030 al (/c./230)2 
0.817 a? (/d:/275)2 

0.719 al 

0.247 al 

0.019 al 
0.022 al (/d/200)2 
0.597 al (/i5j240)2 

0.727 al 

0.042 al 
0.072 al (/r,./230)2 
2.097 al (/i3*/275)2 

1.037 al 



0.30 
0.02 
0.04 
1.03 
0.75 
0.04 
0.03 
0.95 
0.84 
0.29 
0.02 
0.03 
0.70 
0.85 
0.05 
0.08 
2.45 
1.21 



0.31 ±0.04 ±0.02 



0.74 ±0.22 ±0.18 
0.84 ±0.16 ±0.07 



1.14 ±0.42 ±0.28 
0.28 ±0.04 ±0.01 



0.94 ±0.24 ±0.23 
0.73 ±0.15 ±0.03 



2.00 ±0.54 ±0.49 
1.27 ±0.30 ±0.05 



Class II 



i?*V(2S) 



2.262 al 

1.051 al 

3.645 al 

1.939 al 
0.164 al {foh^Qf 
0.230 al {fD'/2S0)^ 
0.111 al (/d/200)2 
0.240 al (/i5./230)2 



0.10 

0.05 

0.16 

0.09 

0.007 

0.010 

0.005 

0.011 



0.075 ±0.021 

< 0.08 
0.153 ±0.028 
0.151 ±0.091 

< 0.033 

< 0.055 

< 0.055 

< 0.117 
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Table 7: Branching ratios (in %) of non-leptonic B decays in the NRSX model 



B- Modes 


NRSX Model af = 1.08 

af = 0.21 


Experimental 
Average 


Class I 






D^D- 

D°D*- 

D^D*- 

D*°D- 

D*°D- 

D*°D*- 

D*°D*- 


0.033 a? (/d/200)2 0.04 
0.938 af(/D, /240)2 1.09 
0.032 a?(/i,./230)2 0.04 
0.873 a? (/z3./275)2 1.02 
0.023 al (/j3/200)2 0.03 
0.639 a?(/D3 /240)2 0.75 
0.077 af (/d./230)2 0.09 
2.235 (/r,|/275)2 2.61 


1.36 ±0.28 ±0.33 
0.94 ±0.31 ±0.23 
1.18 ±0.36 ±0.29 
2.70 ±0.81 ±0.66 


Class II 






K-ip{2S) 

K*-J/iIj 


2.411 al 0.11 
1.122 al 0.05 
3.886 o| 0.17 
2.070 al 0.09 


0.102 ±0.014 
0.070 ±0.024 
0.174 ±0.047 
< 0.30 


Class III 






D*°g- 


0.274 [ai + 1.127 02 (/73/200)]^ 0.48 
0.686 [ai +0.587 02 (fo /200)]^ 0.99 
0.264 [fli + 1.361 a2 {fo- /2:i0)f 0.49 
0.775 [af + 0.661 ai(/D. /230)2 1.19 
+ 1.518 aia2 (/i3-/230)] 


0.50 ±0.05 ±0.02 
1.37 ±0.18 ±0.05 
0.52 ±0.08 ±0.02 
1.51 ±0.30 ±0.02 
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Table 8: Branching ratios (in %) of class I non-leptonic decays in the new model 

described in Section M 



B° Modes 


New Model 


af = 0.98 
af = 0.29 


Experimental 
Average 


Class I 








0.318 aj 


0.30 


0.31 ±0.04 ±0.02 


D+K- 


0.025 al 
0.037 al (/r,/200)2 
1.004 al (/d,/240)2 

0.778 al 


0.02 




D+D- 


0.03 




D+D- 


0.96 


0.74 ±0.22 ±0.18 


D+g- 


0.75 


0.84 ±0.16 ±0.07 


D+K*- 


0.041 al 


0.04 




D+D*- 


0.032 af(/D* /230)2 
0.830 al (/i5./275)2 
0.844 al 


0.03 




D+D*- 


0.80 


1.14 ±0.42 ±0.28 


D+a^ 


0.81 




D*+TT- 


0.296 al 


0.28 


0.28 ±0.04 ±0.01 


D*+K- 


0.022 al 


0.02 




D*+D- 


0.023 al (/i3/200)2 


0.02 




D*+D- 


0.603 al {fDj240f 


0.58 


0.94 ±0.24 ±0.23 


D*+Q- 


0.870 al 


0.84 


0.73 ±0.15 ±0.03 


D*+K*- 


0.049 


0.05 




D*+D*- 


0.085 (/£,*/230)2 
2.414 al (/d;/275)2 
1.217 al 


0.08 




D*+D*- 


2.32 


2.00 ±0.54 ±0.49 




1.16 


1.27 ±0.30 ±0.05 




50.0 aflKfcP 








176.9a2|y„,|2 







30 



Table 9: Branching ratios (in %) of class II non-leptonic -B" decays in the new model. We 
taJie 8 = 20° for the 77-^' mixing angle 



B° Modes 


New Model af = 0.98 
af = 0.29 


Experimental 
Average 


Class II 






LUD*° 


0.800 al 0.07 
0.326 aj 0.03 
2.518 ai 0.21 
1.424 ai 0.12 
0.018 ai 0.002 
0.050 ai 0.004 
0.084 ai(/D/200)2 0.007 
0.116 ai (/d./230)2 0.010 
0.078 ai(/D/200)2 0.007 
0.199 ai(/r../230)2 0.017 
0.081 ai(/D/200)2 0.007 
0.203 ai(/D*/230)2 0.017 
0.058 ai(/D/200)2 0.005 
0.073 ai (/r,./230)2 0.006 


0.075 ±0.021 

< 0.08 
0.153 ±0.028 
0.151 ±0.091 

< 0.006 

< 0.025 

< 0.033 

< 0.055 

< 0.055 

< 0.117 

< 0.057 

< 0.120 

< 0.033 

< 0.050 
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Table 10; Branching ratios (in %) of non-leptonic B decays in the new model 



B- Modes 



New Model 



af = 0.98 
af = 0.29 



Experimental 
Average 



Class I 



D^D- 

D^D*- 
D^D*- 
D*°D- 
D*°D- 
D*'^D*- 
D*°D*- 



0.039 al (/n/200)2 
1.069 al (/r,^/240)2 
0.034 al (/z)./230)2 
0.883 al (/i,./275)2 
0.025 al (/d/200)2 
0.642 al (/d./240)2 
0.091 al (/r,./230)2 
2.570 a? (/zj./275)2 



0.04 
1.03 
0.03 
0.85 
0.02 
0.62 
0.09 
2.47 



1.36 ±0.28 ±0.33 
0.94 ±0.31 ±0.23 
1.18 ±0.36 ±0.29 
2.70 ±0.81 ±0.66 



Class II 



K*-J/ip 
7r~ J/^ 



0.852 al 
0.347 ai 
2.680 al 
1.516 al 
0.038 al 
0.107 al 



0.07 0.102 ±0.014 

0.03 0.070 ±0.024 

0.23 0.174 ±0.047 

0.13 < 0.30 

0.003 0.0057 ±0.0026 

0.009 < 0.077 



Class III 



0.338 [ai +0.729 a2 (/r./200)]^ 
0.828 [ai + 0.450 02 (/d/200)]^ 
0.898 [ai + 0.317 aa (/d/200)]^ 
0.315 [ai + 0.886 aa (/d*/230)]^ 
0.926 [al + 0.456 al [fo- /230)2 

± 1.292 aiaa (/d*/230)] 
1.296 [al + 0.128 (/i3./230)2 

+ 0.269 aia2 (/r../230)] 



£>07r- 
£>*0a]- 



0.48 
1.02 
1.03 
0.48 
1.26 

1.36 



0.50 ±0.05 ±0.02 
1.37 ±0.18 ±0.05 

0.52 ±0.08 ±0.02 
1.51 ±0.30 ±0.02 

1.89 ±0.53 ±0.08 
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ratios are i?i = 1.11 ± 0.23 and i?2 = 1.15 ± 0.34, in agreement with our 
predictions, although the errors are stiU large. We may also consider the 
corresponding ratios for pairs of class 1 transitions to final states which differ 
only in their light meson. For instance, we find 



i?3 = ^-7^ ; f ~ 1.03 [1.06] 



which in the heavy-quark limit become equal unity, too. The experimental 
values are R3 — 0.95 ± 0.24 and R4 — 0.99 ± 0.25, again in agreement with our 
predictions. 

Similar ratios can be taken for class II amplitudes; however, since the 
relevant matrix elements involve current matrix elements between a heavy 
and a light meson, the theoretical predictions are considerably more model 
dependent. Based on the results of our two models, we expect 

i^s = « 0.62 [0.32], 

i?a = - 0.54 [0.23]. (59) 

The strong model dependence does not allow for a test of the factorization 
hypothesis in this case. The corresponding experimental values are R5 = 
0.58±0.11 and Rq — 0.44±0.31, where we have averaged the results for i?° and 
B~ decays (accounting for the different lifetimes) when available. Although 
the experimental errors are sizable, the data seem to prefer the NRSX model. 

We now turn to the actual values of the phenomenological parameters a^^ 
and a|^. From the theoretical point of view, the cleanest determination of af^ 
is from the class of decays B^ D^*^^h~ , where h~ is a light meson {h = tt, 
g or ai). The relevant B^ £)(*)+ transition form factors are known from the 
analysis of semileptonic B decays using the heavy-quark effective theory, and 
the decay constants of the light mesons are experimentally known with good 
accuracy. Also, these transitions have a similar decay kinematics, so that we 
may expect that they are characterized by the same value of . Performing 
a fit to the experimental data, we obtairf[| 

[£,,, = 1.08 ±0.04 [0.98 ±0.04]. (60) 



"^Since we are neglecting final-state interactions, the parameters a'r^ are real numbers, 



and by convention we take a'f^ to be positive. 



33 



The coefficient al^ can also be determined from the decays B 
which are characterized by a quite different decay fcinematics. In principle, it 
would be interesting to investigate whether the resulting value is different in 
the two cases, i.e. whether there is an observable process dependence of the 
phenomenological parameter. In practise, this cannot be done because of the 
large uncertainty in the values of the decay constants of charm mesons. From 
a fit to the data, we find 



.off 



DD,= 1-10 ± 0.07 ± 0.17 [1.05 ± 0.07 ± 0.16] , 



(61) 



where the second error accounts for the uncertainty in /^(*) . In both cases, 

( |60| ) and (|6l|, the data support the theoretical expectation that af^ is close to 
unity [see (p5|)]. 



A value for the parameter 



,cff| 



(but not the relative sign between 



and af^) can be obtained from the class II decays B K^*\j/ip und B 
i4:(*V(2S). From a fit to the six measured branching ratios, we extract 



la. 



eff I 

2 iKii 



= 0.21 ±0.01 [0.29 ±0.01] 



(62) 



This result is more strongly dependent on the form-factor model used, which 
is not surprising given that class II decays involve heavy-to-light transition 
matrix elements. As we have seen above, the NRSX model may be the more 
trustable one; still, we believe the difference between the two results provides 
a realistic estimate of the theoretical uncertainty. 

A determination of from decays with a rather different kinematics is 
possible by considering the class III transitions B~ with h = n or 

g. Moreover, because of the interference of ai and 02 amplitudes, these decays 
are sensitive to the relative sign of the QCD coefficients. From the theoretical 
point of view, it is of advantage to normalize the branching ratios to those of 
the corresponding B^ decays, which are class I transitions. The theoretical 
predictions for these ratios are of the form 



B{B- L>(*)"/i-) _ t(B-) 
B{B° Di*)+h-) ^ r(SO) 



off ^2 eff 



(63) 



where Xi and X2 are process-dependent parameters depending on the ratio of 
some hadronic form factors and decay constants [xi = X2 except for the decay 



B' 



D 



*o , 



For the ratios of branching fractions on the left-hand side, 
we use recent CLEO data reported in Ref. 86, which are more accurate than 
the previous world averages presented in Tables |6|^^10|. Performing a fit to the 
data, we extract the ratio a'^ /a"^ for each channel. The results are collected 
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Table 11: Ratios of non-leptonic decay rates of charged and neutral B mesons and the 
corresponding values for a^^^/aj^ 



Experimental Ratios 



Predictions for Xi 



af/af 



B{B- 



B(BO - 


D+TT-) 


B(S- - 


D^g-) 


B(SO - 




B(S- - 




B(BO - 


D*+TT-) 


B(B- ' 


D*'^g-) 



B{B° D*+g- 



1.73 ±0.25 



= 1.19 ±0.24 



1.64 ±0.28 



1.71 ±0.36 



1.127 [0.729] 

0.587 [0.450] 

1.361 [0.886] 

xi = 0.759 [0.646] 
X2 = 0.813 [0.675] 



0.24 ±0.08 ±0.02 
[0.38 ±0.13 ±0.03] 

0.10 ±0.18 ±0.03 
[0.13 ±0.24 ±0.04] 

0.18 ±0.08 ±0.02 
[0.27 ±0.12 ±0.03] 

0.35 ±0.17 ±0.03 
[0.41 ±0.20 ±0.04] 



in Tahle |Tl|, where the second error results from the uncertainty in the lifetime 
ratioEll t{B-)/t{B°) = 1.06 ± 0.04. Taking the average, and using we 
find 



,eff 



-floh = 0.23 ±0.05 [0.30 ±0.1 



= 0.21 ±0.05 [0.31 ±0.1 



Dh 



(64) 



The value of a|*^ is in good agreement with that obtained in (|6^). Given 
the different decay kinematics in the two processes, this observation is quite 
remarkable. 

The magnitude and, in particular, the positive sign of are of great 
importance for the theoretical interpretation of our results. We find that in 
non-leptonic B decays the two parameters af^ and have the same sign, 
meaning that the corresponding amplitudes interfere constructively. This find- 
ing is in stark contrast tpJLhfi situation encountered in Z3-meson decays, where 
a similar analysis yields 

af [charm = 1-10 ± 0.05 , flf [charm = "0.49 ± 0.04 , (65) 

indicating a strong destructive interference. Since most D decays are (quasi) 
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two-body transitions, this effect is reaponsibie for the observed lifetime dif- 
ference between D+ and D° mesons E3: t{D+)/t{D°) = 2.55 ± 0.04. In B 
decays, on the other hand, the majority of transitions proceeds into multi- 
body final states, and moreover there are many B~ decays (such involving two 
charm quarks in the final state) where no interference can occur. The relevant 
scale for multi-body decay modes may be significantly lower than mt, leading 
to destructive interference (see Figure Therefore, the observed construc- 
tive interference in the two-body modes is not in conflict with the fact that 
t{B-) > t{B°). 

The values for extracted from B — > i^T^*^?/; and B D^*^h decays in 



( |62| ) and (64) indicate that non-factorizable contributions (at the scale fj, — mt,) 
are small in these processes. Using Oj^ji^-^ — C2{mb) + CKipci{mb) = 0.21±0.05 
and a2^\Dh = C2{mb) + CDhCi{mb) = 0.23 ± 0.05 with conservative errors, and 
combining these with the values of the Wilson coefficients given in Table |, we 
find 

Ck4, = 0.44 ± 0.05 , sf^''''^ {rub) 0.11 ± 0.05 , 

Coh = 0.46 ± 0.05 , ei^"'''\mb) = 0.13 ± 0.05 . (66) 

Hence, within errors there is no experimental evidence for a process dependence 
of the value of (, in accordance with our expectation stated in 



9 Tests of Factorization and Extraction of Decay Constants 

Based on the factorization hypothesis, we have made in the previous section 
a number of predictions for the ratios of hadronic decay rates. Within the 
current experimental uncertainties, these predictions agree well with the avail- 
able data. In this section, we shall discuss a particularly, clean method to test 
the factorization hypothesis. As suggested by BjorkenE3, we make use of the 
close relationship between semileptonic and factorized hadronic amplitudes by 
dividing the non-leptonic decay rates by the corresponding differential semilep- 
tonic decay rates evaluated at the same . This method provides a direct test 
of the factorization hypotheais;_moreover, it may be used to determine some 
interesting decay constantsE23ll2J. Assuming that factorization holds, we have 

B(50^^<.).M-) _.e,VI,|.;«l"|K„l'Al;>. (67) 



dB(BO ^ D(*'>+£-v)/dq^' 



with af^ w 1. Here Jm is the decay constant of the meson M, and Vij is the 
appropriate CKM matrix element (depending on the flavour quantum num- 
bers of the meson M). To determine this ratio experimentally, one needs the 
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Table 12: Values of the differential semileptonic branching ratio used in the factorization 



tests, as compiled in Ref. 81 


2 

q 


dB{B 




ml 




0.237 ±0.026 


2 




0.250 ±0.030 


ml. 




0.335 ±0.033 






0.483 ±0.033 


ml. 




0.507 ±0.035 



values of the differential semileptonic branching ratio at various values of q^. 
They have been determined for B D*tv decays in Ref. 81, using a fit to 
experimental data. We collected their results in Table 

Neglecting the lepton mass, we obtain for a pseudoscalar meson P: 



Xp = 



(to| - ml,) 



[ml — {mu + mp)^] [m^ — [mpi — mpY] 



Foiml) 



Fi{ml) 



[mo' + mpY 
1%) 



{mo. - mpY 



m% Ei=o,± \Ht{m%) 



(68) 



The helicity amplitudes Ho{q^) and H±{q^) are defined in the Appendix. For 
the special case that the pseudoscalar meson is a pion, it is an excellent approfc 
imation to expand the quantities and X* in powers of m^/ml, yieldingEZI 



Xt, 



Xl 



Amlmpmn 



l.OOf , 



AmlmBm,p)t 
{m% 



Ami 



m|j,)2 (mB — mD-Y 



0.994. 



Neglecting the tiny deviation from unity, we obtain 

= K = 6^72 l^cff |2 |-^^^|2 ^ |^cff|2 ^ o_gg ggy2 



(69) 



(70) 



This prediction may be compared with the experimental value i?* = 1.18 ± 
0.22 GeV^ obtained by combining the — > D*~^tt~ branching ratio from 
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Table ^ with the value for the differential semileptonic branching ratio at = 
given in Table This yields = 1.11 ± 0.10, in good agreement with 
the expectation based on factorization. An even cleaner test of factorization is 
obtained when ( |67| ) ia_^valuated for a vector or pseudovector meson, in which 
case one has exactly E 
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Xv = X*y^l. (71) 

Since the lepton pair created by the {V — A) current carries spin one, its 
production is kinematically equivalent to that of a (pseudo-) vector particle 
with four- momentum g^. For a g meson in the final state, for instance, we 
thus obtain 

Re = Rl^ 6^V' l«f P \Vud? ~ \af\' X 2.48 GeV^ , (72) 

to be compared with the experimental value R* — 2.92 ± 0.71 GeV^. This 
gives al^ = 1.09 ± 0.13, again in good agreement with the expectation based 
on factorization. In principle, eqs. ( [70| ) and ( |7^ ) offer the possibility for four in- 
dependent determinations of the QCD parameter a\^. Good agreement among 
the extracted values supports the validity of the factorization approximation 
in B decays. Already at the present level of accuracy, it shows that there is 
little room for final-state interactions affecting the magnitude of the considered 
decay amplitudes. 

From the kinematic argument about the equivalence of the lepton pair in 
the semileptonic decay and the spin-1 meson in the hadronic decay, it follows 
that ( |67| ) is valid separately for longitudinal and transverse polarization of the 
D* meson in the final state. Thus, the polarization of the D* meson produced 
in the non-leptonic decay D*~^V~ should be equal to the polarization 

in the corresponding semileptonic decay B D*£D at — my. However, 
in order to turn this prediction into a test of the factorization hypothesis one 
would have to determine the polarization of the D* meson with high precision. 
This is so because in the semileptonic as well as in the non-leptonic case the D* 
polarization at the points q'^ ~ and q^ ~ Qmax determined by kinematics 
alone to be 100% longitudinal and 1/3 longitudinal, respectively. This shows 
that for a stringent test of the factorization hypothesis at small one must 
determine the transverse polarization contribution with a small relative error. 
Especially for the semileptonic decays, a precision measurement of the q^ de- 
pendence of the polarization appears to be a complicated task, however. Still, 
we can make a rather precise prediction for this quantity using heavy-quark 
symmetry. In the heavy-quark limit, the ratio of transverse to longitudinal 
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Table 13: Theoretical predictions for the ratio T'Y/Ttot at fixed g^, where Ftot = Tt + Fl 





ml ml^ ml, q^^^^ 


rx/rtot 


12 ±1 26 ± 2 48 ±1 2/3 



polarization at some fixed is simply given by 

Ft Aq^ {m?g + m\,, — q^ 



Fl {mB ~ mD*Y[{mB + mo'Y - q^] ^'^'^^ 
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Including the leading symmetry-breaking corrections to this result 
obtains the numbers shown in Table For the polarizatiorijof the D* meson 
in the decay B° D*+ g- , the CLEO Collaboration finds Fx/Ftot = (7 ± 
5 ± 5)%, in agreement with our prediction of 12% transverse polarization for 
the semileptonic decay at q^ ~ ml. However, in order for this test to be 
sensitive to deviations from factorization, the experimental uncertainty will 
have to be reduced substantially . The situation may be more favourable 
in the decay with predicted 48% of transverse polarization, 

hopefully allowing for a measurement with smaller relative uncertainties. 

We shall now discuss an alternative use of (|67|). Assuming the validity 
of the factorization hypothesis with a fixed value for a^;^, one mn.y.. ertiploy 



this relation for the determination of unknown decay constants F H i. In 



particular, we can use it to determine the decay constants of the Ds and D* 
mesons. As explained above, in the latter case we simply have Xd* = X'^, — 1. 
But also in the case of mesons we obtain essentially model-independent 
predictions. Defining the mass ratio x — ^{mj^i,) + m£i^)/mB, we find in the 
heavy-quark limit Ell (x w 0.36 and x w 0.38, respectively) 

l-3x^ + 2x^ ^ 



1 — T — 9t^ 

Xh. - " ^ 0.37. (74) 

1 — x + 2x^ 

These values are close to the predictions of the NRSX model, which are 

«1.33, X*o^k0.39. (75) 

Using these values (allowing for a theoretical error of ±0.03) together with 
|a°^| = 1.08 ± 0.04, we obtain the theoretical predictions: 

Rd, = (87.2 ±6.7)/!,^ , 
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- (25.6 ±2.7)/!,^, 
Rd^ - = (65.6 ± 4.8) /I,. . (76) 

By averaging the experimental data on and B~ decays into two charm 
mesons (taking into account the differences in the i3-meson hfetimes), we ob- 
tain 

B{B^DD-) = (0.95 ±0.24)%, 
B{B~>DD*-) = (1.00 ±0.30)%, 
B{B^D*D-) = (1.03 ±0.27)%, 

B{B ^ D*D*-) = (2.26 ± 0.60)% . (77) 

Using the last two averaged decay rates gives the experimental ratios 

R*o^ = (2.13 ± 0.58) GeV^ , i?!,. = (4.46 ± 1.22) GeV^ . (78) 

Comparing this with the theoretical predictions in ([76|), we find 

fo, = (288 ± 42) MeV , /^j = (261 ± 36) MeV . (79) 

More precise determinations of the decay constants are possible if, instead 
of using (|6^), we consider ratios of non-leptonic decay rates, comparing pro- 
cesses involving Ds and D* mesons with those involving the light mesons tt and 
Q. These processes involve a similar kinematics, so that the ratios of the corre- 
sponding decay rates are sensitive to the same form factors, however evaluated 
at different values. This method has the advantage that in the ratios the 
phenomcnological parameter cancels; similarly, we may hope that some of 
the experimental systematic errors cancel. Using the NRSX model, we find 



B(B0 


- D+D-) 


B(BO 


D+n-) 


B{B^ - 




B(BO - 


D*+iT-) 


B(^o - 




B(B" - 


-^D+Q-) 


B(^o ^ 





1.01 



0.72 



= 0.74 



B(B" ^ D*+g-) V /< 



1.68 I ^ ) . (80) 



I e 



40 



Combining these predictions with the average experimental branching ratios 
in (77), we find the rather accurate values 



fa, = (234 ± 25) MeV , /d; = (271 ± 33) MeV . (81) 

The result for fo^ is in excellent agreement with the value fo^ = 241 ±37 MeV 
in (^), extracted from leptonic decays of Ds mesons. The ratio fo'/fD^ — 
1.16 ± 0.19, which cannot be determined,iDm leptonic decays, is in good 
agreement with theoretical expectations l 
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Along these lines, there are numerous other possibilities for extracting 
information on decay constants and current matrix elements. In particular, 
the decay constants of P-wave particles like oq, ai, Kq and Ki are of interest 
due to their sensitivity to relativistic quark motion, which allows for a test of 
hadron models. For instance, one may use the ratio 

BiBO^D*-a^)^^_^^(f^y (82) 



B(i?o ^ D*+g-) • V /, 



to determine the pseudovector meson decay constant fa^ , which is defined by 
( I Qlp-lsQ l^n) — ^p.'m-ai fai ■ From the experimental results, we obtain 

fa, = (1.22 ± 0.19) fe = (256 ± 40) MeV , (83) 

which agrees with the large value derived from r decays shown in Table ^ 
The fact that the decay constant of a P-wave meson is of the same order of 
magnitude as those of the corresponding S-wave mesons is quite remarkable. 
It shows that mesons containing light quarks are highly relativistic systems; 
for non-relativistic constituent quarks one would expect the wave function of 
a state with non-zero orbital momentum to vanish at the origin, leading to a 
vanishing decay constant. 



10 B Decays to Baryons 

We end our discussion of non-leptonic processes with a remark on the par- 
ticularly interesting case of B decays into baryon-antibaryon pairs. For these 
transitions it is probably not sufficient to apply factorization in the way de- 
scribed so far. The relevant flavour flow diagram (for a 6 ^ c transition) is 
shown in Figure ^ In contrast with the situation realized in B decays into 
mesons, the c and d quarks produced in the weak interaction now end up in 
the same hadron. It is, therefore, appropriate to rewrite the effective weak 
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Figure 3: Flavour flow diagram for a B decay into a baryon-antibaryon pair. 



Hamiltooia 



ordering til 



in a different form using charge- conjugate fields and a Fierz re- 



c^i^i) {cd')liubh^+c+if,) (cd'iiub), 
penguin operators , 



h.c. 



(84) 



wliere (cd')3 



£kij {I — 75) is an {S — P) colour-antitriplet diquark 



current. The colour-sextet current {cd')Q is given by a similar expression. The 
appearance of the coefficients c^{^j,) and c+(/i) shows that perturbativc QCD 
interactions affect the colour-antitriplet and sextet channels differently: the 
QCD force is attractive (c_ > 1) in the antitriplet and repulsive {c^ < 1) in 
the sextet channel. 

From the very existence of baryons it is evident that, besides the pertur- 
bativc enhancement, there must also be a strong long-distance binding force 
between two quarks in a colour-antitriplet state. This can be taken into ac- 
count by introducing diquark fields and their couplings to the local currents 
occurxing in (84). QCD sum rules predict that these diquark couplings are 
large Ei3. This result has been confirmed by the jinagnitude of AS* = 1 transi- 
tions, where diquark effects play a dominant roletd. Thus, we expect B decays 
to baryon-antibaryon pairs to proceed predominantly via the formation of a 
diquark-antidiquark state, followed by the creation of a quark-antiquark pair 
in the colour field of these diquarks. This decay mech anism i mplies interesting 
selection rules, which allow for a number of predictionsllij'tij. For instance, the 
charm baryon in baryonic 6 — > c transitions is built from (cd) spin-zero states. 
Thus, in B decays one should not find a (cuu) nor a (cms) nor a (ess) baryon, 
as long as final-state interactions may be neglected. In 6 — > u transitions one 
naturally obtains a A/ = i selection rule, because a spin-zero (ud) diquark 
necessarily has isospin zero. No A resonances should therefore be produced in 
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these decays, but there is no such suppression of A production. A comparison 
of these predictions with experimental data, once these will become available, 
should provide information about the creation process of quark-antiquark pairs 
inside hadrons. For an estimate of the expected branching ratios, the reader 
is referred to Ref. 114. 

11 Summary 

We have presented an overview of the theory and phenomenology of exclu- 
sive hadronic decays of B mesons, concentrating on two-body modes. Such 
decays are strongly influenced by the long-range QCD colour forces. Theoret- 
ically, their description involves hadronic matrix elements of local four-quark 
operators, which are notoriously difficult to calculate. The factorization ap- 
proximation is used to relate these matrix elements to products of current 
matrix elements. Conventionally, the factorized decay amplitudes depend on 
two phenomenological parameters ai and 02, which are connected with the 
Wilson coefficients Ci{fi) appearing in the effective weak Hamiltonian. We 
have shown that (except for decays into two vector mesons) this approach can 
be generalized in a natural way to include the dominant non-factorizable con- 
tributions to the decay amplitudes. In the generalized factorization scheme, 
the effective parameters af^ and become process-dependent. However, 
using the large- TVc counting rules of QCD we have argued that in energetic 
two-body decays of B mesons af^ ?a 1 and al*^ ~ C2{mb) + (ci{mb), where 
C = 0{1/Nc) is a dynamical parameter. Moreover, we have shown that the 
process dependence of ( is likely to be very mild, so that it can be taken to be 
a constant for a wide class of two-body decays. These theoretical expectations 
are supported by the data. From a fit to the world average branching ratios of 
two-body decay modes, we obtain af^ « 1.08 and « 0.21, corresponding to 
C pa 0.45. There is no evidence for a process dependence of these parameters; 
in particular, the values obtained for from the decays B — » K'^*'^'ij) and 
B- D(*)0/i-, where h = n OT g, are in good agreement with each other. 

The obvious interpretation of the fact that the ratio /af^ is positive and 
the value of the parameter ( close to the "naive" factorization prediction ( = 
1/3 is that, in energetic two-body decays of B mesons, a fast moving colour- 
singlet quark pair interacts little with soft gluons. Hence, factorization works 
at a high scale of order mf,. The situation is different from that encountered 
in the much less energetic decays of D mesons, where one finds a negative 
value of /af^, corresponding to C ~ 0. Since the energy release in charm 
decays is much less, even soft gluons can rearrange the quarks, and the effective 
factorization scale is lower. According to the relation between the factorization 
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scale and the ratio of the phenomenological parameters exhibited in Figure ^ 
this leads to a negative value of a"^ /a'^ and thus a smaller value of Q. 

The most important ingredient of factorized decay amplitudes are the 
hadronic form factors parametrizing the hadronic matrix elements of quark 
currents. In the case of the heavy-to-heavy transitions B D and B — > 
D* , heavy-quark symmetry implies simple relations between the various form 
factors. Incorporating the leading symmetry-breaking corrections using the 
heavy-quark effective theory, it has become possible to extract all B — > Z?*^*^ 
form factors from semileptonic decay data with good precision. The use of 
heavy-quark symmetry constitutes a significant improvement over earlier esti- 
mates of non-leptonic amplitudes, which were based on model calculations of 
the relevant form factors. As a consequence, for all class I decays considered 
in this article the factorized decay amplitudes can be predicted without any 
model assumptions. 

There has not been similar progress in the calculation of current matrix 
elements between heavy and light mesons. For these transitions we must still 
rely on phenomenological models. Consequently, the theoretical predictions 
for class II decay amplitudes involve larger theoretical uncertainties. In order 
to get an idea about the amount of model dependence, we have considered two 
different quark models and compared their results. We find that it is possible 
to determine the parameter with a theoretical accuracy of about 25%. 
In the context of each model, a large set of class II branching ratios can be 
reproduced within the experimental errors using a fixed value of . 

We have discussed various tests of the (generalized) factorization hypoth- 
esis by considering ratios of decay rates, and by comparing non-leptonic decay 
rates with semileptonic rates evaluated at the same value of q^. Within the 
present experimental uncertainties, there are no indications for any deviations 
from the factorization scheme in which a^^ and are treated as process- 
independent hadronic parameters. Accepting that this scheme provides a use- 
ful phenomenological concept, exclusive two-body decays of B mesons offer 
a unique opportunity to measure the decay constants of some light or charm 
mesons, such as the ai, Ds and D*. In particular, on the basis of the ex- 
perimental data available today, we find /^^ = (234 ± 25) MeV and fo* — 
(271 ± 33) MeV. The result for is in excellent agreement with the value 
extracted from leptonic decays of Ds mesons. The ratio fn* / fo^ = 1.16 ±0.19 
is in good agreement with theoretical expectations. 
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Appendix 

We collect our definitions for the weak decay form factors, which parametrize 
the hadronic matrix elements matrix elements of flavour-changing vector and 
axial currents between meson states. For the transition between two pseu- 
doscalar mesons, Pi{p) P2{p'), we define 

{P,{p')\ V,\P,{p)) = ({p + p'), - ^^^i^ g,) F,{q') + I^i-^ q,Fo{q') , 

(A.l) 

where = {p — p')^ is the momentum transfer. For the transition of a pseu- 
doscalar into a vector meson, Pi{p) V2{e,p'), we define 

{V2{e,p')\V^\Mi{p)) = e^.^ffe*VVT^(9^), 
nil + 



{V2{e,p')\A^\Mi{p)) 



(mi + ma) e* Ai{q^) ^— ^ {p + p')^^ ^2(9^) 

^ 1711+1712 



21712 ^^qj.Asiq'^) 



+ 2m2^-^qy,Ao{q^), 



(A.2) 



where is the polarization vector, satisfying e • p' = 0. Here, the form factor 
^43(5^) is given by the linear combination 

Moreover, in order for the poles at = to cancel, we must impose the 
conditions 

Fi(0) = i^o(0), A3(0) = Ao(0). (A.4) 
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The helicity amplitudea_£|o('j'^) and H±{q^) are given by the foUowmg combi- 
nations of form factors tij {K denotes the momentum of the daughter meson 
in the parent rest frame): 



H±{q^) = (mi +m2) =p 
1 



miK 



nil + '7i2 
{■m\ - ml - q^){mi + Ai{q^) 



mi + 7712 



(A.5) 
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